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ABSTRACT 

The  pulsed-gradient  spin-echo  (PGSE)  NMR  technique  was  used  to  acquire  a  broad  range  of 
solvent  self-diffusion  coefficients  for  the  toluene/polystyrene  (MW  "570  000)  system.  Self 
diffusion  coefficients  were  determined  for  samples  having  polymer  weight-fractions  ranging  from 
0.04  to  0.90  over  a  temperature  range  of  25-115  °C.  Several  different  approaches  to  interpretation 
of  the  diffusion  data  have  been  considered.  The  temperature  dependence  of  the  data  exhibits 
Arrhenius  behavior  with  energies  of  activation  that  increase  with  concentration  from  2.6  to 
16.2  kcal/mole.  Several  theories  based  on  the  assumption  that  the  solvent  diffusion  is  slowed  by 
the  polymer  obstructing  its  diffusion  path  are  also  used  to  interpret  this  data.  Our  results  support 
the  hypothesis  that  the  normalized  diffusion  coefficient  is  relatively  independent  of  thermodynamic 
parameters.  The  data  are  in  qualitative  agreement  with  a  "standard  curve"  for  normalized  diffusion 
coefficients.  This  approach  shows  potential  as  a  tool  for  predicting  diffusion  rates  over  a  wide 
concentration  range  for  a  variety  of  polymer  solvent  systems.  The  final  approach  used  to  interpret 
this  data  is  that  of  the  free  volume  theory  proposed  by  Vrentas  and  Duda.  The  agreement  between 
the  free  volume  theory  and  experimental  data  is  improved  considerably  if  two  of  the  needed 
parameters  are  adjusted  to  give  optimum  fits  to  the  experimental  diffusion  data.  This  approach 
yielded  a  root  mean  squared  of  the  ^ctional  difference  between  experiment  and  theory  of  14%. 
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Introduction 

Small  molecule  transpon  in  polymer  solutions  is  of  critical  imponance  in  many  processes.  It 
can  determine  polymerization,  curing  and  drug  delivery  rates;i  plasticization,  pigment  stability, 
diy'ing  and  coloration  of  coatings;2  and  many  other  properties  of  polymeric  systems.  The  role  of 
diffusion  in  most  of  these  processes  has  been  reviewed  elsewhere.3-6  Although  diffusion  in 
gasses  and  small  molecule  liquids  is  better  understood,  there  remains  considerable  uncertainty 
about  such  processes  in  polymer  solutions.  Much  of  this  uncertainty  is  a  result  of  the  variety  of 
different  approaches  for  the  determination  of  diffusion  coefficients.  Different  approaches  yield 
different  diffusion  coefficients  whose  inter-relationship  is  unclear.  The  different  diffusion 
coefficients  and  the  various  methods  for  determining  them  in  polymer  systems  have  also  been 
reviewed.2  Self-diffusion  is  a  measure  of  transport  due  to  Brownian  motions  of  molecules  in  the 
absence  of  a  chemical  concentration  gradient.  One  could  argue  that  the  presence  of  the  nuclei, 
labelled  by  their  precessional  frequencies  (as  in  the  PGSE  experiment),  creates  a  thermodynamic 
concentration  gradient.  However,  these  gradients  do  not  significantly  affect  the  Brownian  motion 
in  the  sample  and  the  measurement  is  relatively  non-perturbing.  The  lack  of  a  chemical 
concentration  gradient  simplifies  development  of  a  theoretical  understanding  of  this  transport 
process.  Thus  self-diffusion  can  serve  as  a  starting  point  for  the  analysis  of  the  diffusion 
processes  in  polymeric  systems.  In  the  past,  seif-diffusion  measurements  were  made  using  radio¬ 
isotope  tracers.  Such  procedures  were  long  and  tedious  because  of  the  synthetic  pathways.  As  a 
result,  the  amount  of  self-diffusion  data  for  polymer  systems  found  in  the  literature  is  moderate  at 
best.  With  the  advent  of  the  pulsed-gradient  spin-echo  (PGSE)  NMR  technique,  self-diffusion 
coefficients  can  be  routinely  determined  in  favorable  circumstances  in  about  thirty  minutes.  It  is 
therefore  feasible  to  perform  a  complete  study  of  solvent  self-diffusion  in  an  ideal  system. 

In  the  present  report,  we  have  chosen  to  study  toluene  diffusion  in  solutions  of  polystyrene 
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because  this  system  has  been  studied  in  great  detail  and  many  of  the  parameters  needed  in  the 
analysis  of  the  diffusion  data  have  already  been  determined.  Our  analysis  will  begin  with  the  more 
simple,  and  empirical  approaches  to  interpretation  of  self-diffusion  data.  Arrhenius  plots  of 
diffusion  coefficients  will  be  used  to  determine  activation  energies  of  the  diffusion  process  as  a 
function  of  concentration.  We  will  then  consider  an  analysis  based  on  obstruction  theories  and 
modified  Enskog  theory  which  corrects  for  momentum  transfer  between  the  solvent  and  polymer. 
Finally,  the  free  volume  theory  of  Vrentas  and  Duda  (V&D)  will  be  considered.  The  advantages 
and  disadvantages  of  each  approach  will  be  discussed. 
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EXPERIMENTAL 

The  PGSE  technique  used  was  similar  to  the  Stilbs  modification*  of  the  Stejskal-Tanner^ 
method.  The  technique  uses  a  pulsed  magnetic  field  gradient  to  create  phase  variarions  of  nuclear 
spins  which  are  determined  by  the  position  of  the  nuclei  in  the  field.  Diffusion  results  in 
incomplete  rephasing  of  the  spins  following  the  application  of  a  refocussing  pulse.  The  signal 
intensity  is  exponentially  related  to  the  diffusion  coefficient  and  the  length  of  the  gradient  pulse  in 
an  isotropic  sample.  By  taking  several  spectra  with  different  gradient  pulse  lengths,  one  can 
extract  a  diffusion  coefficient  from  the  rate  of  signal  decay  as  a  function  of  the  gradient  pulse 
length.  Reviews  of  the  techniquei0-i2  and  its  application  to  polymeric  systems'*. have  been 
published.  One  limitation  of  the  technique  is  that  the  signal  intensity  is  attenuated  by  spin-spin 
relaxation  ( T2  )  as  well  as  the  diffusion  coefficient.  At  ambient  temperatures,  samples  with  very 
high  polymer  concentration  (above  80%  w/w  polymer)  had  T2's  which  were  prohibitively  short. 
Thus  it  was  not  possible  to  determine  self-diffusion  coefficients  in  all  samples  at  all  of  the 
temperatures  of  interest.  Also,  since  the  signal  to  noise  ratio  degrades  as  the  relaxation  time 
becomes  shorter,  the  precision  of  the  data  for  the  high  concentration  samples  was  not  as  good  as 
that  of  the  other  samples. 

The  diffusion  coefficients  were  determined  from  the  slope  of  a  linear  least-squares  fit  of  the 
log  of  peak  height  vs  P  =  62(A-  5/3);  where  5  is  the  gradient  pulse  width  and  A  is  the  time 
between  the  initial  pulse  and  the  refocussing  pulse.'O  For  each  diffusion  measurement,  eleven 
spectra  were  typically  obtained  using  different  gradient  pulse  widths.  Often,  a  few  of  the  spectra 
taken  at  long  gradient  widths  were  discarded  because  of  insufficient  signal  to  noise.  At  the  value 
of  A  used  in  this  study  (typically  67  ms),  the  contributions  of  the  polymer  protons  to  the  toluene 
methyl  resonance  is  negligible  because  of  the  short  T2  of  the  polymer  signal.  The  standard 
deviation  of  the  slope  of  the  log-linear  least  squares  fit  was  taken  as  an  indication  of  the  precision 
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of  the  diffusion  measurement  and  ranged  from  1  to  6%  of  the  measured  diffusion  coefficient. 
However,  the  accuracy  of  the  measurement  in  the  PGSE  experiment  was  estimated  to  be  about 
10%  of  the  diffusion  coefficient  being  measured  in  the  range  studied.  The  field  gradient  strength 
(0.047  T/m)  was  determined  using  neat  cyclohexane  as  a  reference  material.  The  value  of  the 
cyclohexane  self-diffusion  coefficient  was  taken  to  be  1.47  x  10-9  mffs  at  25°C. 

All  experiments  were  performed  on  a  JEOL  FX-90Q  spectrometer  operating  at  89.5  MHz  for 
protons.  The  spectra  were  taken  using  an  external  lithium  lock  so  it  was  not  necessary  to 
incorporate  deuterated  solvent  into  the  samples.  All  spectra  were  taken  using  a  1000  Hz  window, 
4000  points  and  a  pulse  delay  of  8  seconds.  The  number  of  scans  required  to  attain  reasonable 
signal  to  noise  ranged  from  8  to  128.  A  JEOL  variable  temperature  control  unit  was  used  in  the 
elevated  temperature  studies  and  was  capable  of  maintaining  the  temperature  to  ±  2  °C.  The 
temperature  controller  was  calibrated  with  a  blank  sample  prior  to  the  diffusion  measurements. 

Samples  were  prepared  by  placing  weighed  amounts  of  the  constituents  in  5  mm  NMR  tubes. 
The  tubes  were  then  sealed  and  placed  in  a  heating  block  at  ca.  60  °C  until  homogeneous. 
Equilibration  times  were  on  the  order  of  several  weeks  for  the  most  concentrated  samples.  Toluene 
and  polystyrene  were  purchased  from  Aldrich  Chemical  Co.,  Milwaukee,  WI  and  Scientific 
Polymer  Products  Inc.,  Ontario,  NY,  respectively.  Both  were  used  as  received.  The  polymer 
molecular  weight  was  characterized  by  gel  permeation  chromatography  (GPC)  and  was  found  to 
have  an  average  molecular  weight  of  270,000  g/mol  and  a  polydispersity  of  1.5.  It  has  been 
demonstrated  that  the  molecular  weight  distrubution  does  not  affect  the  solvent  self-diffusion 
coefficient  in  polymer-solvent  systemsj^  Therefore  this  broad  molecular  weight  distribution  is 
acceptable  for  the  present  study. 

Polymer  volume  fractions  were  calculated  from  weight  fractions  using  the  tabulated 
temperature  dependencies  of  the  densities  of  neat  toluene'5  and  polystyrenei^  assuming  ideal 
volumes  of  mixing.  The  densities  vary  from  0.862  to  0.773  g/ml  and  from  1.079  to  1.029  g/ml 
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over  the  temperature  range  studied  for  toluene  and  polystyrene,  respectively. 
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Results  and  Discussion 

Toluene  self-diffusion  coefficients  were  determined  over  the  temperature  range  of  25  to  115  °C 
at  ten  degree  intervals  for  fifteen  samples  having  concentrations  ranging  from  3  to  89  %  w/w 
polymer.  The  results  of  these  measurements  are  presented  in  Table  I  as  a  function  of  weight 
fraction  polymer.  The  lack  of  data  for  the  concentrated  polymci  samples  at  low  temperature  was 
due  to  the  limitations  imposed  by  short  spin-spin  relaxation  times  (T2)  due  to  the  relatively  slow 
molecular  motion  in  those  samples.  A  sudden  decrease  in  T2  at  very  high  concentrations  and  low 
temperatures  was  observed  and  probably  was  due  to  the  proximity  of  the  glass  transition  (Tg)  in 
these  samples.  The  concentration  dependence  of  Tg  for  the  toluene/PS  system  has  been  previously 
studied  for  a  system  of  MW  =  18,000  g/mol.!"^  Interpolation  of  this  data  indicates  that  the  Tg  of 
samples  with  polymer  weight  fractions  of  0.791, 0.830,  0.852  and  0.886  should  be  -1,  7,  17,  and 
35  °C  respectively.  17  However,  the  molecular  weight  of  the  polymer  used  in  that  analysis  was 
considerably  below  the  value  at  which  the  limiting  Tg  is  expected.  Therefore,  we  can  assume  that 
in  the  samples  used  in  the  present  study,  the  Tg  will  be  somewhat  higher  than  that  predicted  in  the 
above  analysis.  For  each  of  these  samples  the  lowest  temperature  at  which  the  diffusion  coefficient 
could  be  determined  was  approximately  50  ®C  above  the  Tg.  This  observation  suggests  a 
relationship  between  free  volume  which  determines  the  polymer  molecular  motions  and  the  solvent 
diffusion.  Also,  since  a  change  in  the  dynamics  of  the  system  occurs  at  the  glass  transition 
temperature,  diffusion  measurements  below  this  point  would  be  indicative  of  different  processes 
than  those  observed  in  solution.  Thus  the  data  set  presented  here  is  a  complete  representation  of 
solvent  self-diffusion  in  the  liquid  like  solurion  over  the  given  temperature  range. 

Self-diffusion  coefficients  for  neat  toluene  were  also  determined  and  listed  in  Table  I.  The 
diffusion  coefficients  and  energy  of  activation  (see  Table  II)  for  diffusion  of  neat  toluene  are 
consistent  with  data  reported  in  the  literature  (Ej  =  2.5  kcal/mol).i8.i9  The  overall  trends  in  the 
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data  for  the  polymer  solutions  were  consistent  with  one’s  expectations;  diffusion  coefficients 
increase  as  temperature  increases  and  decrease  as  polymer  concentration  increases.  The  solution 
data  were  also  quantitatively  consistent  with  the  limited  self-diffusion  data  that  has  been  previously 
reported20-23  for  this  system.  Several  approaches  to  interpreting  this  data  are  available  bui  only  a 
few  will  be  considered  for  simplicity,  insight  gained  and  predictive  ability. 

Arrhenius  Behavior 

Diffusion  coefficient  data  are  often  analyzed  in  terms  of  an  Arrhenius  energy  of  activation.  In 
dilute  solutions,  where  the  diffusing  species  does  not  often  encounter  a  polymer  molecule,  the 
diffusion  rate  is  limited  by  the  energy  required  for  the  diffusing  species  to  escape  its  present 
surroundings  and  move  into  an  adjacent  environment.  Thus,  in  this  regime  one  would  expect  the 
Arrhenius  approach  to  be  successful.  However,  in  moderate  to  high  concentration  solutions  the 
diffusion  process  is  limited  by  the  polymer  molecular  motions.  Such  motions  are  highly 
cooperative,  involving  several  monomer  units  within  the  chain  as  well  as  those  polymer  segments 
and  solvent  molecules  adjacent  to  the  chain  undergoing  the  motion.  One  would  not  necessarily 
expect  such  a  process  to  be  characterized  by  a  single  energy  of  activation.  However,  there  has 
been  some  experimental  evidence  that  supports  the  validity  of  the  Arrhenius  type  of  approach. 
Mutual-diffusion  coefficients  of  toluene  in  polystyrene  were  shown  to  exhibit  Arrhenius  behavior 
at  weight  fractions  of  0.10  and  0.04  and  over  the  temperature  range  of  1 10-178  Although 
self  and  mutual-diffusion  are  formally  different  processes,25  they  involve  similar  mechanisms. 
Self-diffusion  data  have  also  been  presented  which  show  Arrhenius  behavior  for  this  and  similar 
systems  over  limited  temperature  range  at  high  polymer  concentrations.22 

The  logarithm  of  the  toluene  self-diffusion  coefficient  vs  inverse  temperature  is  plotted  in 
Figure  1.  The  weight  fraction  of  polymer  in  each  sample  is  listed  in  the  legend.  The  deviations 
from  linearity  apparent  in  the  data  arc  attributed  mainly  to  inaccuracies  of  the  temperature 
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controller.  The  larger  amount  of  scatter  in  the  low  concentration  samples  may  have  been  due  to  the 
uncontrolled  effects  of  convection.  In  the  more  concentrated  samples  convection  did  not  play  a 
significant  role  because  of  the  higher  viscosity.  The  data  in  Figure  1  are  linear  within  experimental 
error. 

The  energies  of  activation  (Ea)  and  intercepts  (  D°  =  D  at  infinite  temperature)  are  listed  in 
Table  II  and  were  determined  from  log  linear  least-squares  fit  of  the  plots  in  Figure  1.  The 
standard  deviations  of  the  linear  fits  are  also  listed  in  the  table.  Activation  energies  increase 
gradually  with  concentration  in  the  dilute  regime  and  increase  more  rapidly  at  concentrations  above 
60  %  w/w  polymer.  The  values  of  D°  show  a  similar  trend.  A  plot  of  energy  of  activation  as  a 
function  of  concentration  is  shown  in  Figure  2.  The  error  bars  represent  the  standard  deviations  in 
the  linear  fits  used  to  calculate  tl.  ^  activation  energies.  Energies  of  activation  for  toluene  diffusion 
in  polystyrene  have  previously  been  determined  for  samples  having  concentrations  in  the  range  of 
80-95%  polymer  for  the  temperature  range  of  60-105  °C.22  These  data  have  been  included  in 
Tabie  II  and  Figure  2  for  comparison  and  are  similar  to  those  in  this  work.  A  gradual  increase  in 
activation  energy  followed  by  a  marked  increase  at  high  polymer  concentration  is  apparent. 
However,  the  transition  to  rapidly  changing  activation  energy  occurs  at  slightly  lower  polymer 
concentration  in  the  present  work.  In  addition,  these  data  are  also  consistent  with  lower 
concentration  diffusion  data  for  benzene  and  cyclohexane  in  polystyrene  solutions. 

The  concentration  dependence  of  the  glass  transition  temperature  (Tg)  in  polystyrene  solutions 
of  molecular  weight  similar  to  that  used  in  the  previous  study22  has  been  analyzed  for  a  variety  of 
solvents  including  toluene.  12  The  Tg  for  a  95%  PS/toluene  system  was  reponed  to  be  63  °C.  This 
corresponds  with  the  highest  concentration  and  lowest  temperature  at  which  diffusion 
measurements  were  reported  for  the  18,000  PS  system.  NMR  diffusion  measurements  can  not  be 
determined  below  the  Tg  on  our  system  because  of  the  very  shon  spin-spin  relaxation  times  in 
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those  samples.  Our  most  concentrated  samples  were  glasses  at  room  temperature  as  indicated  by 
stress  cracking  of  the  samples  on  cooling  to  room  temperature. 

The  differences  in  the  behavior  of  the  energies  of  activation  in  our  two  data  sets  is  likely  due  to 
the  lower  molecular  weight  of  polymer  used  in  the  previous  studies.22  The  relationship  between 
glass-transition  temperature  and  molecular  weight  has  been  previously  reported  for  polystyrene. 
This  study  indicated  that  the  lower  molecular-weight  bulk  polymer  has  a  glass  transition 
temperature  that  is  significantly  lower  (87  °C)  than  the  limiting  value  (100  °C)  expected  for  the 
polymer  used  in  the  present  study.  The  apparent  relationship  between  activation  energy  for  solvent 
diffusion  and  Tg  suggests  that  the  same  free  volume  may  determine  both.  The  lower  Tg  of  the  low 
molecular  weight  polymer  could  be  interpreted  as  an  indication  of  a  larger  amount  of  free  volume 
available  in  the  bulk  material.  The  activation  energy  is  relatively  independent  of  concentration  until 
the  available  free  volume  reaches  a  critical  level.  At  concentrations  above  this  level,  stronger 
temperature  dependencies  are  observed.  The  concentration  at  which  the  transition  occurs  in  the 
low  molecular  weight  polymer  solutions  occurs  at  a  higher  concentration  because  the  low 
molecular  weight  polymer  has  more  free  volume  than  the  higher  molecular  weight  polymer. 

These  observations  suggest  that  temperatures  should  be  referenced  to  the  Tg  when  performing 
calculations  of  activation  energies.  This  approach  has  been  applied  with  some  success  to  mutual- 
diffusion  data  at  zero  solvent  fraction  using  free  volume  theories.^f'  Extension  of  that  analysis  to 
self-diffusion  is  beyond  the  scope  of  this  report.  In  the  less  concentrated  solutions,  which  are  well 
above  Tg,  the  diffusion  coefficient  measurements  were  independent  of  the  polystyrene  molecular 
weight.23 

Although  the  Arrhenius  analysis  shows  limited  predictive  capability,  it  does  offer  some  insight 
into  the  mechanism  of  diffusion.  In  solutions  of  low  to  moderate  concentration  (less  than  50% 
polymer)  the  energy  of  activation  is  only  slightly  dependent  on  concentration.  This  implies  that  the 
transport  process  is  limited  by  a  single  step  with  a  measurable  activation  energy.  In  contrast,  the 
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rapidly  changing  activation  energy  observed  in  the  concentrated  samples  implies  a  more 
complicated  process.  These  observations  are  consistent  with  free  volume  theories  which  predict 
that  in  dilute  solution  the  diffusion  coefficient  is  determined  by  the  interactions  of  the  solvent  with 
its  environment.  However,  in  concentrated  solutions  the  transport  is  determined  by  the  availability 
of  free  volume.  The  movement  of  free  volume  is  a  cooperative  process  and  is  not  expected  to  be 
characterized  by  a  single  activation  energy. 

Microscopic  Approaches 

There  are  several  common  approaches  for  treating  the  microscopic  behavior  of  solvent 
transport  in  polymer  solutions.  Several  of  these  have  been  reviewed  for  diffusion  in  gels  by  Muhr 
and  Blanshard.27  Formally,  a  loosely  crosslinked  gel  and  polymer  solution  are  similar  with 
respect  to  solvent  transport  because  the  local  segmental  motions  of  the  polymer  to  which  the 
solvent  transport  would  be  sensitive  are  similar.  This  has  previously  been  shown  for  swollen 
polystyrene  beads.28  It  is  the  microscopic  solvent  environment  which  controls  solvent  motion 
rather  than  the  more  macroscopic  motion  of  entire  polymer  chains  which  determine  the 
macroscopic  viscosity.  Muhr  and  Blanshard  cite^^:  i)  the  obstruction  effect;  ii)  increased 
hydrodynamic  drag;  iii)  alteration  of  solvent  properties  (e.g.  intermolecular  spacings);  and  iv) 
polymer  involvement  as  major  contributions  for  reduced  solvent  diffusion  coefficients  as  compared 
to  bulk  solvent. 

Since  the  polymer  moves  much  slower  than  the  solvent  it  acts  as  an  obstruction.  The  solvent 
must  diffuse  around  the  relatively  stationary  polymer,  thereby  increasing  its  diffusion  path  and 
decreasing  its  diffusion  coefficient.  The  degree  of  obstruction  is  determined  by  the  volume  fraction 
of  the  polymer  in  the  solution.  Obstruction  theories  have  been  very  popular  because  of  their 
simplicity.  For  the  present  work,  we  are  content  to  discuss  only  two  of  these  namely  the  Maxwell- 
Fricke  (or  simply  the  Maxwell  model)  and  Mackie-Mears  models. 
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The  Maxwell27.  29, 30  model  for  diffusion  is  formally  similar  to  that  for  electric  conductivity 
and  permeability.  It  assumes  immobile,  impenetrable  species  suspended  in  a  mobile  solvent 
continuum.  It  has  been  shown  that  the  Maxwell  theory  predicts  the  maximum  reasonable 
diffusion  coefficient  in  the  presence  of  obstructions  when  the  obstructions  are  assumed  to  be 
spheres.27  This  approach  has  been  shown  to  accurately  predict  the  diffusion  coefficient  of  water  in 
suspensions  of  latex  and  silica  particles  at  low  concentrations.^!.  32  This  is  taken  as  an  indication 
that  the  theory  accounts  for  the  effects  of  obstructions  on  the  diffusion  rate  for  large  particles.  If 
one  assumes  that  this  theory  holds  even  for  relatively  small  polymer  segments  which  are  stationary 
relative  to  the  solvent,  the  Maxwell  equation  for  diffusion  simplifies  to:27. 3i,  32 

D(l-<l>)/Do  =  (l-<t>')/(l-h<pVx)  (1) 

where  D  is  the  observed  diffusion  coefficient,  Dq  is  the  diffusion  coefficient  of  the  neat  solvent,  (}) 
is  the  volume  fraction  of  polymer,  ([)'  is  the  volume  fraction  of  the  polymer  plus  any  "non¬ 
diffusing"  solvent  and  %  is  the  shape  factor.  Non-diffusing  solveiit  is  defined  as  that  just  next  to  a 
polymer  segment  with  a  very  slow  diffusion  coefficient.  The  shape  factor  can  take  values  ranging 
from  2  for  spheres  to  1.5  for  rods.  We  will  assume  a  value  of  2  for  X-  For  the  discussion  below, 
the  difference  between  1.5  and  2  for  x  is  small  compared  to  the  differences  between  their 
predictions  and  the  data.  For  this  model  the  (1  -  <j) )  correction  term  has  been  applied  to  DA^q  in 
order  to  account  for  the  proper  reference  volume  in  the  solution  to  Pick's  law.3i-34  if  as  a  first 
approximation  we  set  <()  =  (()',  equation  (1)  reduces  to: 

D/Do  =  1/(1 +(1)/X)  (2) 

This  has  the  physical  significance  of  assuming  that  there  is  no  non-diffusing  solvent  which  is 
appropriate  for  colloidal  species  with  minimal  surface  effects.^!.  32 

A  second  approach  by  Mackie  and  Mears27. 35  js  based  on  tortuosity  and  a  simple  cubic  lattice 
model  with  similar  size  sites  for  the  solvent  and  polymer  segments  (repeat  units).  A  more  general 
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theory  based  simply  on  the  volume  fraction  of  polymer,  <j)  (rather  than  the  fraction  of  repeat  units) 
yields: 

D/Do  =  (1-<|))2/(1 +0)2  (3) 

This  model  deals  with  the  small  size  of  the  impenetrable  polymer  segment.  Extensions  of  this 
model  have  been  suggested  which  consider  the  redistribution  of  polymer  segments.  For  rapid 
polymer  segment  redistribution,  this  results  in27 

D/Do  =  (1-0)2  (4) 

Slower  molecular  motion  should  yield  a  result  between  these  two  extremes.  This  approach  is 
interesting  in  light  of  studies  which  suggest,  through  a  simple  free  volume  treatment,  that  the 
segmental  motion  of  the  polymer  and  diffusion  of  the  solvent  are  correlated.  This  has  been  shown 
in  the  poly(isopropyl  acrylate) — chloroform  system^S  but  is  believed  to  be  generalizable  to  other 
polymer  solutions. 

The  modified  Enskog  result  is  based  on  kinetic  theory  and  accounts  for  momentum  transfer 
between  rigid  spheres.37  The  polymer  is  modelled  as  a  set  of  immobile  beads.  The  details  of  this 
approach  will  be  discussed  elsewhere.  With  the  modified  Enskog  theory,  the  solvent  self¬ 
diffusion  ceofficient  is  given  by:^^ 

1/D=  l/D,,p-t- 1/Dj5  (5) 

where  Dj  p  and  D^  ^  are  the  pair  (mutual)  diffusion  coefficients  of  the  tracer  relative  to  the  polymer 
and  other  solvent  molecules,  respectively.  The  value  of  the  Dj's  are  calculated  from  equation  6  in 
reference  37.  The  ratio  of  the  molecular  diameters  of  the  solvent  to  that  of  the  polymer  (mer) 
segment  is  a  fitted  parameter. 

Shown  in  Figure  3  are  the  self-diffusion  coefficients,  normalized  by  Dp,  as  a  function  of 
polymer  volume  fraction.  Also  shown  in  the  figure  are  the  predictions  of  the  Maxwell,  Mackie  and 
Meares  and  modified  Enskog  theories.  At  first  glance  the  most  striking  feature  of  Figure  3  is  that 
the  normalized  diffusion  coefficients  fall  on  very  similar  curves.  In  fact,  solvent  self-diffusion 
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coefficients  from  a  number  of  systems  also  fall  within  experimental  error  of  these 
curves.23. 33. 38, 39  A  Standard  plot  of  this  type  can  be  used  to  predict  diffusion  rates  in  a  variety 
of  systems  with  a  knowledge  of  the  densities  of  the  constituents  and  the  self-diffusion  coefficient 
of  the  neat  solvent.  This  type  of  behavior  is  anticipated  when:38  the  system  is  well  above  Tg 
(T  >  Tg  50°C);  the  polymer  molecular  weight  is  high;  and  in  the  case  of  multi-solvent  systems, 
when  there  are  no  strong  specific  polymer-small  molecule  interactions.  Lee  and  Lodge  have  found 
that  methyl  red  in  the  methyl  red-toluene-poly(vinyl  acetate)  system  does  not  follow  this  universal 
behavior  because  it  preferentially  binds  to  poly(vinyl  acetate)  and  toluene  does  not.38 

A  definite  temperature  dependence  is  found  in  the  normalized  diffusion  coefficients  in  Figure  3. 
In  general,  the  normalized  diffusion  coefficient  (D/Dq)  increases  with  increasing  temperature, 
although  these  effects  are  not  much  larger  than  the  experimental  errors.  Similar  effects  are  also 
found  in  this  system  at  lower  temperatures,  even  if  volumes  of  mixing  are  accurately  accounted 
for.40  The  magnitude  of  this  effect  is  small  which  is  not  surprising  because  the  energies  of 
activation  are  similar  from  0-50  wt  %  polymer.  Unfortunately,  with  thiS  type  of  analysis  all  of  the 
data  depend  on  the  accuracy  of  the  Dq  measurement  which  is  difficult  to  make  because  of 
convection. 

A  comparison  of  the  theories  with  the  experimental  values  in  Figure  3  suggest  that  the  Maxwell 
theory  gives  a  poor  approximation  to  the  experimental  data,  while  the  others  seem  more 
reasonable.  The  experimental  data  fall  well  below  the  Maxwell  theory.  This  could  be  due  to 
setting  (])  =  ({>'.  An  alternative  approach  using  the  Maxwell  theory  would  be  to  use  equation  1  and 
solve  for  (])'  for  a  given  value  of  D(1  -  <!))/Do.  The  difference  between  <})  and  <J)'  could  then  be  used 
to  estimate  the  amount  of  non-diffusing  solvent. 

In  aqueous  solutions  it  is  common  to  assume  that  the  non-diffusing  solvent  arises  because  of 
the  association  of  water  molecules  with  the  polymer  through  hydrogen  bonding.  This  approach 
has  been  used  to  measure  hydration  numbers  in  solutions  of  water  soluble  polymers33  and  in 
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biological  systems.34,41,42  in  many  cases  the  limiting  form  of  the  Maxwell  theory  is  used  at  low 
concentrations.  A  comparison  of  solvent  self-diffusion  coefficients  in  a  variety  of  polymer-solvent 
systems  has  shown  that  the  apparent  volume  of  the  polymer  bound  solvent  is  relatively 
independent  of  the  presence  of  hydrogen  bonding,  solvent  quality,  polymer  molecular  weight  and 
temperature.23.38  This  implies  that  the  chemical  interaction  between  the  solvent  and  polymer  is  not 
a  critical  factor  in  determining  the  diffusion  coefficient.  Others  have  noted  that  the  normalization  of 
the  diffusion  coefficient  as  is  generally  required  by  the  obstruction  theories  removes  many 
thermodynamic  effects  from  the  diffusion  data.20.43 

We  believe  that  the  interaction  between  the  diffusing  species  and  polymer  which  causes  the 
disagreement  between  the  Maxwell  theories  and  experimental  data  is  physical  (not  chemical)  in 
nature.  By  this  we  mean  that  specific  chemical  interactions  (e.  g.  hydrogen  bonding)  do  not 
themselves  cause  the  difference.  We  have  previously  called  this  the  "surface  drag"  effect  which 
has  been  briefly  discussed. 27  We  envision  that  the  solvent  molecules  next  to  the  polymer  molecule 
move  at  roughly  the  same  rate  as  the  polymer.  Since  the  polymer  segments  move  very  slowly  the 
solvent  next  to  them  also  moves  very  slowly  compared  to  the  bulk  solvent.  There  may  in  fact  be 
several  layers  of  solvent  molecules  which  are  motionally  restrained  by  the  polymer.  Dielectric 
dispersion  data  for  solutions  of  polyethylene  oxide  in  water  at  microwave  frequencies  show  a 
broad  range  of  water  relaxation  times,  ranging  from  one  to  ten  times  that  of  pure  water.44  This 
indicates  that  water  is  present  in  several  different  states  in  contrast  to  the  two  state  model  assumed 
in  the  calculation  of  hydration  numbers.  The  two  state  model  does,  however,  have  the  advantage 
of  being  computationally  and  conceptually  simple.  The  chemical  interaction  between  the  solvent 
and  polymer  does  determine  the  lifetime  of  the  solvent  next  to  a  polymer  segment;  however,  this 
time  is  very  short  relative  to  the  diffusion  time  in  the  PGSE  experiment  even  in  strongly  associated 
systems.  The  solvent  exchanges  rapidly  between  the  different  sites  during  the  experiment  and 
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therefore  one  observes  an  average  diffusion  coefficient.  The  diffusion  rate  is  determined  only  by 
the  number  of  sites  for  the  solvent  interaction  and  the  diffusion  rate  of  the  neat  solvent.  This 
behavior  is  not  observed  in  suspensions  of  latex  paiticles^i.  32  because  the  surface  area  of  particles 
is  very  small  relative  to  the  number  of  sites  in  the  polymer  solution. 

The  prediction  of  the  Mackie  and  Meares  theory  is  much  closer  to  the  experimental  data  than  the 
Maxwell  theory  but  falls  slightly  below  it.  The  limiting  form  of  this  model,  assuming  fast  polymer 
motion,  brackets  the  data  on  the  other  side  (not  shown).  It  seems  reasonable  tha.  intermediate 
molecular  motion  would  fit  the  data  somewhere  in  between  these  two  extremes.  At  present  there  is 
no  simple  extension  of  the  theory  to  accommodate  some  limited  polymer  motion,  but  the  effect  of 
temperature  could  be  explained  by  the  decrease  in  the  relative  rates  of  solvent  to  polymer  motion. 

The  modified  Enskog  theory  fits  the  data  quite  well.  This  model  properly  accounts  for  the 
momentum  transfer  between  the  solvent  and  polymer.^s  The  theory  was  fit  to  the  data  with  one 
parameter,  the  radius  of  a  polymer  segment,  or  equivalently,  the  ratio  of  the  solvent  to  polymer 
segment  diameter.  The  curve  shown  is  fit  with  this  ratio  equal  to  0.9.  We  expect  this  treatment  to 
not  be  valid  at  high  polymer  concentrations,  ie.  >  50  %.  This  type  of  analysis  can  predict  the 
temperature  dependence  as  well  and  we  are  currently  extending  this  treatment  to  other  systems.'^s 

These  types  of  analyses  are  useful  in  the  study,  of  a  wide  variety  of  polymer  solvent  systems 
and  over  a  wide  range  of  temperatures  and  concentrations  below  ca.  60  %  polymer.  However,  at 
concentrations  above  approximately  60%  w/w  polymer  the  magnitude  of  the  experimental  errors  is 
larger  than  that  of  the  normalized  diffusion  coefficient.  Unfortunately,  it  is  in  this  concentration 
regime  that  many  polymer  processes  are  diffusion  limited.  Thus  an  alternative  approach  to  the 
analysis  of  this  data  is  considered  in  the  next  section. 

Free  Volume 


Free  volume  may  be  defined  as  the  volume  of  a  system  which  is  not  accounted  for  by  the 


Pickup  and  Blum,  Self-Diffusion  of  Toluene  in  Polystyrene  Solutions 


18 


volume  of  the  molecules  as  determined  by  the  van  der  Waals  radii  of  the  atoms  that  make  up  the 
molecules.  It  is  assumed  that  this  additional  volume  is  distributed  randomly  in  the  system  and  that 
it  rearranges  with  little  or  no  activation  energy.  When  a  portion  of  this  free  volume  of  sufficient 
size  moves  to  a  location  adjacent  to  a  diffusing  molecule,  the  molecule  may  jump  into  that  space. 
The  diffusion  process  is  assumed  to  be  a  series  of  these  jumps. 

Cohen  and  TumbulH^  first  proposed  the  dependence  of  diffusion  coefficients  on  the  available 
free  volume  for  the  case  of  hard  spheres.  This  theory  was  adapted  for  polymer  solutions  by 
Fujita'*'^.  48  and  was  shown  to  be  reasonably  successful  at  predicting  self-diffusion  coefficients 
over  limited  concentration  and  temperature  ranges.36. 49. 50  This  theory  was  funher  modified  by 
Vrentas  and  Duda  (V&D)  in  order  to  account  for  the  differences  in  the  volumes  required  for 
diffusion  by  the  solvent  and  polymer  jumping  units.51  It  has  been  shown  that  the  Fujita  theory  is  a 
special  case  of  the  more  general  V&D  theory  in  which  the  solvent  and  polymer  jumping  units 
require  the  same  free  volume  to  jump.52  Experimental  evidence  has  been  presented  .vhich 
correlates  the  molar  volume  of  the  jumping  unit  of  a  diffusing  solvent  with  the  diffusion 
coefficient.23. 53  xhe  limited  success  of  the  Fujita  theory  is  attributed  to  the  fact  that  the  systems 
studied  meet  the  criterion  of  similar  size  for  the  solvent  and  polymer  jumping  units.  In  addition, 
the  dependence  of  the  Fujita  theory  on  the  reference  diffusion  values  magnifies  the  importance  of 
the  accuracy  of  the  chosen  reference.  When  the  extrapolated  solvent-diffusion  coefficient  is  used 
as  the  reference,  this  problem  may  be  even  worse.54 

Another  free  volume  theory  proposed  by  Paul55  is  similar  to  that  of  V&D  in  its  dependence  on 
free  volume,  however  it  includes  a  factor  to  account  for  the  possibility  of  the  solvent  jumping  back 
to  the  position  from  which  it  originated  but  ignores  the  effects  of  the  relative  ^izc  of  solvent  and 
polymer  jumping  units.  The  predictions  of  the  Paul  and  V&D  theories  have  been  compared.56 
The  results  suggested  that  the  Paul  theory  is  reasonably  accurate  at  low  to  moderate  concentrations 
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(where  the  jump-back  factor  is  expected  to  play  a  significant  role)  but  is  ineffective  at  high  polymer 
concentrations.  Consequently  we  have  chosen  the  free  volume  theory  proposed  by  V&D  for 
comparisons  to  our  experimental  results. 

The  V&D  theory  states  that  the  self-diffusion  coefficient  is  given  by;‘'i 

D  =D*  e.xp{-Ea/RT}  exp{-Y  (coiVi*  )A^fh}  (6) 

where  D*  is  a  fitted  parameter,  Ea  is  the  activation  energy  for  moving  the  diffusing  molecule  from 
its  initial  environment  to  its  destination  (not  the  same  energy  of  activation  discussed  in  section  "a"), 
R  is  the  gas  constant,  T  is  temperature,  y  accounts  for  the  fact  that  the  same  free  volume  is  available 
to  more  than  one  molecule  for  diffusion  and  may  take  values  ranging  from  1/2  to  1,  V*  is  critical 
volume  (the  volume  of  the  diffusing  species)  of  the  solvent  or  polymer  at  0  K,  Vfh  is  the  average 
free  volume  available  for  diffusion  in  the  system,  ^  is  the  ratio  of  the  solvent  critical  free  volume  to 
that  of  the  polymer  jumping  unit  and  accounts  for  the  fact  that  the  monomer  is  not  the  actual 
jumping  unit,  co  is  the  mass  fraction,  and  the  subscripts  1  and  2  refer  to  the  solvent  and  polymer, 
respectively.  The  ratio  Y  A^fh  is  can  be  detemiined  from  available  free  volume  parameters  in  the 
following  manner,56 

VFH/Y=(Kii/Y)a)i(K2i+T-Tgi)-F  (  K12/Y ) t02 (  K22  +  T  -Tg2)  (7) 

where  T  is  temperature,  Tg  is  the  glass  transition  temperature,  CD  is  the  weight  fraction,  and  the 
subscripts  1  and  2  refer  to  solvent  and  polymer  respectively.  The  K's  are  constants  which  are 
evaluated  as: 


K„/y=2V°i  a, 

(8) 

K12/Y  =  V  V(2.303  Cl  C2) 

(9) 

K22  =  C2 

(10) 

where  V®i  is  the  molar  volume  of  the  solvent  at  the  glass  transition  temperature,  tt]  is  the 
expansion  coefficient  for  the  solvent,  and  the  Cs  are  the  Williams,  Landel  and  Ferry*^  (WLF) 
constants  where  the  subscripts  indicate  which  constant  is  used.  Since  the  value  of  (K22  -  Tg2)  for 
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this  system  is  negative  (Table  HI),  equation  6  suggests  that  whenever  the  temperature  falls  below 
(K22  -  Tg2),  the  polymer  makes  a  negative  contribution  to  the  available  free  volume.  However,  it  is 
unrealistic  to  expect  that  the  addition  of  more  material  would  remove  free  volume.  In  our  analysis, 
this  was  compensated  for  by  setting  the  polymer  contribution  to  free  volume  equal  to  zero 
whenever  the  temperature  was  below  this  critical  value.  This  correction  to  the  theory  was  found  to 
have  a  minimal  etfect  on  the  outcome  of  the  analysis. 

The  theory  does  have  a  large  number  of  parameters;  however,  many  of  them  can  be  determined 
independently.  The  critical  free  volumes  for  diffusion  are  assumed  to  be  equal  to  the  molar 
volumes  at  0  K  and  were  calculated  using  the  methods  discussed  by  Haward.57  Values  of  0.917 
and  0.850  cmVg  were  used  for  the  critical  free  volumes  of  the  solvent  and  polymer  at  0  °K, 
respectively.  The  WLF  constants  are  determined  from  viscosity  data  and  have  been  tabulated  for  a 
variety  of  polymer  solvent  systems.'*^,  50  a  preliminary  analysis  of  the  data  showed  that  the 
contribution  of  the  Ea  term  to  the  quality  of  the  fit  was  minimal.58  Therefore  Ea  was  set  to  zero  and 
the  remaining  two  parameters,  ^  and  D* ,  were  adjusted  to  minimize  the  root  mean  of  squared  of 
the  fractional  difference  between  the  theory  and  experimental  data  which  is  given  by  the  following 
equation: 


%  deviation  =  100  x  /V  ^  X  (11) 

where  Dex  and  Dcak  are  the  experimental  and  theoretical  self-diffusion  coefficients  and  the 
summation  is  over  the  n  experimental  data  points.  This  type  of  analysis  is  preferable  to  standard 
least  squares  fits  because  the  data  cover  two  and  one-half  orders  of  magnitude.  The  values  of  the 
optimized  parameters  as  well  as  the  remaining  parameters  used  in  this  analysis  are  listed  in 
Table  ni.  Plots  of  the  experimental  data  and  the  predictions  of  equation  6  are  shown  in  Figure  4. 
The  solid  lines  in  the  figure  represent  the  theoretical  predictions.  Only  a  portion  of  the 
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experimental  data  is  shown  in  the  figure  in  order  to  minimize  overlap  although  the  entire  data  set 
was  used  in  the  analysis.  The  quality  of  the  fit  of  the  remaining  data  was  similar  to  that  shown  in 
the  figure.  These  values  resulted  in  an  average  difference  squared  between  theory  and  experiment 
of  ca.  14%.  We  note  that  the  present  data  is  in  better  agreement  (lower  absolute  error)  than 
previously  reported  data  from  both  self  and  mutual-diffusion  measarements.56  In  that  study  the 
mutual  diffusion  coefficients  were  converted  to  self-diffusion  coefficients.  The  best  fit  values  of  ^ 
and  D*  found  in  this  study  are  similar  to  those  reported  previously.56  The  value  of  D*  is  about 
20%  different  and  probably  due  to  the  different  range  of  temperatures  studied  in  this  work.  ^  ^ 

The  overall  fit  of  the  experimental  data  and  theory  is  good;  however,  a  systematic  variation 
between  the  two  is  apparent.  Specifically,  in  the  low  concentration  regime  an  upward  curvature  is 
apparent  in  the  experimental  data.  This  curvature  could  not  be  accounted  for  by  the  theory  as 
presented  here.  This  behavior  has  been  observed  in  other  systems  and  has  been  addressed 
previously.59  It  has  been  suggested  that  this  behavior  is  due  to  changes  in  the  value  of  E,  as  a 
function  of  temperature.  In  the  low  concentration  region,  the  environment  of  the  solvent  molecules 
is  changing  rapidly.  The  environment  is  dominated  by  solvent  molecules  and  not  polymer  chains. 

Thus  the  energy  of  interaction  between  the  solvent  and  its  environment  is  expected  to  change.  At 
present  the  details  of  this  interaction  have  not  yet  been  determined  and  their  relationship  to  system 
conditions  is  unclear.  A  further  analysis  of  this  interaction  would  offer  greater  utility  to  the  free 
volume  approach  to  solvent  self-diffusion. 

Finally,  we  should  point  out  that  similar  types  of  measurements  and  comparisons  have  been 
made  using  forced  Rayleigh  scattering^^'^  or  photobleaching  recovery  experiments.  In  principle 
these  measurements  may  be  applicable  over  a  wider  range  than  the  data  presented  here,  but  depend 
on  the  presence  of  a  probe  molecule.  In  polystyrene-  tetrahydrofuran  solutions,  the  self-diffusion 
of  methyl  red^®  was  found  to  be  adequately  fit  by  cither  the  Fujita  or  Vrentas-Duda  theory.  This  is 
in  agreement  with  the  present  work. 


/ 


Pickup  and  Blum,  Self-Diffusion  of  Toluene  in  Polystyrene  Solutions 


22 


Conclusions 

The  usefulness  of  any  theory  for  diffusion  is  determined  by  its  ability  to  predict  experimental 
phenomena  and  the  amount  of  insight  into  the  diffusion  process  that  is  gained  by  its  application. 
Each  of  the  three  approaches  presented  here  is  somewhat  successful  in  this  regard.  When  treated 
in  an  Arrhenius  format,  the  self-diffusion  data  yield  energies  of  activation  as  a  function  of 
concentration.  It  is  apparent  from  such  data  that  the  diffusion  process  changes  dramatically  as  one 
approaches  the  glassy  polymer.  The  behavior  of  activation  energy  at  low  polymer  concentration 
suggests  that  the  diffusion  process  can  be  defined  by  a  single  activation  energy  and  perhaps  a 
single  step.  At  concentrations  above  60%  w/w  polymer,  the  activation  energy  changes  rapidly 
with  concentration  which  implies  a  multi-step  or  cooperative  process.  Though  such  an  analysis 
has  little  predictive  capability,  it  does  yield  some  insight  into  the  diffusion  process.  With  some 
modification,  such  that  activation  energies  are  referenced  to  glass  transition  temperatures,  such  an 
analysis  could  show  improved  predictive  capability. 

Three  different  approaches  to  obstruction  theory  were  considered.  The  simplified  Maxwell 
theory  is  by  far  least  accurate  at  prediction  of  diffusion  coefficients  in  polymer  solutions.  This  is  at 
least  in  part  due  to  the  assumption  that  there  is  no  solvent  associated  with  the  polymer  to  cause  an 
added  obstruction  volume.  However,  the  Maxwell  theory  is  useful  in  that  it  predicts  a  maximum 
possible  diffusion  coefficient  in  these  systems  and  shows  the  obstruction  effect  alone  can  not 
predict  the  experimental  results.  The  second  approach  to  obstruction  theories  to  be  considered  was 
that  of  Mackie  and  Meats  which  is  based  on  a  cubic  lattice  model.  This  approach  gives  reasonable 
agreement  but  falls  slightly  below  the  experimental  data.  Finally,  the  modified  Enskog  theory 
accounts  for  momentum  transfer  and  fits  the  data  within  experimental  error. 

The  obstruction  theories  have  strong  predictive  abilities  since  plots  of  the  normalized  diffusion 
coefficient  are  relatively  independent  of  all  parameters,  including  the  system  being  studied.  Thus  a 
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standard  curve  could  be  used  to  predict  diffusion  rates  with  a  knowledge  of  the  densities  of  the 
components  and  the  diffusion  coefficient  of  the  neat  solvent.  This  data  is  generally  available,  but 
the  approach  is  limited  to  polymer  concentrations  below  about  60%.  The  discrepancy  between 
theory  and  experiment  is  accounted  for  by  a  surface-drag  effect.  Though  there  is  some  evidence  to 
support  the  surface  drag  concept,  it  has  yet  to  be  directly  confirmed. 

Free  volume  theories  describe  diffusion  on  a  microscopic  level  and  give  some  details  about  this 
process.  Though  a  large  number  of  parameters  are  needed  to  fit  the  theory,  these  parameters  are 
available  for  a  variety  of  polymer-solvent  systems.  In  our  analysis  ,two  parameters  were  adjusted 
to  give  best  fits  of  the  experimental  data.  It  is  possible  to  approximate  the  adjustable  parameters 
from  limited  diffusion  data  in  order  to  give  the  theory  predictive  capability.  This  interpretation  has 
strengths  in  both  prediction  and  understanding  the  diffusion  process.  Systematic  differences 
between  the  theory  and  experiment  suggest  that  further  modification  of  the  theory  is  in  order. 
Development  of  a  specific  relationship  between  the  interaction  energy  and  concentration  would 
significantly  improve  the  agreement  between  theory  and  experiment. 
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Table  I.  Toluene  self-diffusion  coefficients  a  in  polystyrene  solutions 
function  of  concentration  and  temperature. 


Temoerature 

Conc.c 

25 

35 

45 

55 

65 

75 

85 

95 

105 

115 

0 

26.6 

29.4d 

32.2 

34.5 

45.4 

48.3 

55.3 

58.7 

61.6 

72.7 

0.036 

24.1 

26.9 

29.7 

30.5 

41.1 

43.4 

46.8 

53.3 

59.3 

66.9 

0.150 

19.0 

21.2 

23.2 

25.7 

30.3 

35.4 

38.7 

44.4 

46.8 

53.1 

0.201 

16.9 

18.8 

20.1 

23.0 

28.1 

33.4 

34.8 

41.9 

44.2 

49.8 

0.267 

13.6 

15.2 

17.4 

19.6 

23.9 

27.8 

31.1 

34.9 

37.7 

41.2 

0.369 

9.5 

11.2 

12.6 

14.2 

18.4 

20.4 

23.5 

26.9 

31.4 

33.0 

0.444 

7.26 

8.27 

9.46 

10.3 

14.3 

16.6 

19.1 

22.1 

24.9 

27.0 

0.486 

5.97 

6.85 

8.15 

8.78 

12.9 

13.7 

16.1 

19.8 

21.7 

24.2 

0.580 

3.15 

3.89 

4.65 

5.54 

7.57 

9.41 

11.0 

13.0 

15.0 

17.7 

0.603 

2.57 

3.15 

4.08 

4.66 

6.58 

8.:  8 

9.55 

11.4 

13.9 

15.6 

0.661 

1.51 

1.85 

2.30 

2.94 

4.50 

5.70 

6.84 

8.18 

10.2 

11.7 

0.714 

0.79 

0.95 

1.37 

1.61 

2.80 

3.43 

4.49 

5.66 

7.04 

8.54 

0.791 

0.51 

0.59 

0.60 

1.06 

1.28 

1.77 

2.55 

3.31 

3.97 

0.830 

0.43 

0.64 

0.87 

1.17 

1.89 

2.49 

0.852 

0.27 

0.44 

0.78 

1.21 

2.08 

0.886 

0.27 

0.46 

0.85 

a  All  diffusion  coefficients  are  in  units  of  10  -10  m2/sec.  b  Temperature  was  maintained 
with  a  precision  of  ±  2  °C.  <=  Polymer  weight  fraction.  ^  Interpolated  from  pure  toluene 
data. 
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Table  II.  Calculated  energies  of  activation  for  diffusion  of  toluene  in 
polystyrene  solutions. 

weight  fraction  PS  (kcal/mole)  D°  (m2/s)  a 


0.0 

2.61  ±0.15 

-15.4  ±0.2 

0.036 

2.66  ±0.15 

-15.4  ±0.4 

0.150 

2.72  ±0.10 

-15.5  ±0.1 

0.201 

2.89  ±0.11 

-15.4  ±  0.2 

0.267 

2.98  ±  0.07 

-15.4  ±0.1 

0.369 

3.33  ±  0.09 

-15.2  ±0.1 

0.444 

3.60  ±0.15 

-15.0  ±  0.2 

0.486 

3.77  ±0.15 

-14.9  ±0.2 

0.580 

4.53  ±  0.10 

-14.2  ±0.1 

0.603 

4.77  ±0.10 

-14.1  ±0.1 

0.661 

5.5210.16 

-13.3  ±0.2 

0.714 

6.40  ±  0.20 

-12.5  ±0.3 

0.791 

6.61  ±  0.42 

-13.1  ±0.6 

0.830 

9.'7±0.33 

-10.2  ±0.5 

0.852 

13.66  ±  0.30 

-4.6  ±  0.4 

0.886 

16.26  ±  0.65 

-2.1  ±0.9 

0.789 

8.3b 

0.850 

8.4b 

0.890 

9.0b 

0.939 

16.5b 

a  Dq  is  reported  in  units  of  10  -10  m2/sec.  bprom  reference  22.  See  text  for  details. 
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Table  III.  Values  for  Free  Volume  Parameters  used  to  Fit  the  Diffusion  Data. 


Kii/Y 

1.57  X  10-3  cm3/gKa 

K12/Y 

5.82  X  10  cm3/g  K  a 

K21  -Tgi 

-90.5  K  a 

K22  -Tg2 

-327  K  a 

Vi* 

0.917  cm3/g 

V2* 

0.850  cm3/g 

Ea 

0 

0.536  b 

D* 

3.85  X  10  cmVsec  b 

a  from  reference  55,  ^  determined  by  optimization  of  die  d.^ta  presented  here. 


Log  D  (niVsec) 


Figure  1.  Plots  of  log  D  vs  1/T  for  toluene  in  polystyrene  at  several  different  weight  fractions 
Oisted  in  the  legend).  The  lines  drawn  are  best  fits  and  indicate  the  energy  of  activation. 
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Figure  2.  Energies  of  activation  for  toluene  seif-diffusion  in  polystyrene  solutions  plotted  as  a 
function  of  concentradon.  The  curves  are  drawn  to  aid  the  eye  only. 


Volume  Fraction  Polymer 


Figure  3.  Normalized  plots  of  toluene  diffusion  coefficients  in  polystyrene  at  several  different 
temperatures  (figure  legend)  as  a  function  of  composition.  Also  shown  are  the  predictions  of  the 
Maxwell  ( _ );  Mackie  and  Meares  (— )  and  modified  Enskog  (•••)  theories. 
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